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Abstract 


We propose a new model to create a light meson in the heavy quarkonium transition, which is 
inspired by the Nambu—Jona-Lasinio (NJL) model. Hadronic transitions of JP? = 177 higher char- 
monia with the emission of an 7 meson are studied in the framework of the proposed model. The model 
shows its potential to reproduce the observed decay widths and make predictions for the unobserved 
channels. We present our predictions for the decay width of Y > J/wn and © > h,(1P)n, where Y 
are higher S and D wave vector charmonia, which provide useful references to search for higher char- 
monia and determine their properties in forthcoming experiments. The predicted branching fraction 
B(w(4415) > h.(1P)n) = 4.62 x 1074 is one order of magnitude smaller than the J/yn channel. Esti- 
mates of partial decay width r(Y — J/wn) are given for Y (4360), Y (4390) and Y (4660) by assuming 
them as cé bound states with quantum numbers 33 D1, 33D: and 5°91, respectively. Our results are in 
favor of these assignments for Y (4360) and Y (4660). The corresponding experimental data for these 
Y states has large statistical errors which do not provide any constraint on the mixing angle if we 
introduce S— D mixing. To identify Y (4390), precise measurements on its hadronic branching fraction 
are required which are eagerly awaited from BESIII. 


1 Introduction 


Quantum chromodynamics (QCD), the gauge theory of strong interaction, received huge devolvement 
during the last few decades. However, it is still a subject of intensive research of various theoretical 
constructs (for instance, see recent review [1]). Study of heavy quarkonium decays is a good probe to 
understand the nonperturbative nature of QCD at different energy regimes. Thanks to the wealth of 
experimental facilities like CLEO, Belle, BABAR, CDF, DO, and BESIII, now we have intensive experi- 
mental data in the charmonium (c2) energy regime. This provides us with great opportunities to test and 
explore the nature of strong interactions in the heavy quark sector. Currently BESIII is taking data in 
the cé energy regime and it is easy to produce JPC = 177 higher charmonia through e+e~ annihilations. 
Figure 1 is the sketch of the intermediate production of vector charmonia, which further decay into J/q7. 
Since the center-of-mass energy of BESIII can go up to 4.6 GeV [2], which is around the mass region 
of (5S) and u(4D), it is a good opportunity to study the production and decay mechanism of higher 
vector charmonia. In the future PANDA also plans to collect data [3] in the cé@ energy regime which 
ete” colliders are not capable of producing directly. These experimental facilities will surely help us to 
deepen our understanding of heavy quarkonium physics and hence the nonperturbative aspects of strong 
interaction. 
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Figure 1: Feynman diagram associated with the intermediate production of higher vector charmonia via 
ete” annihilation, which further decay into J/wn. 


Heavy hadron spectroscopy has celebrated almost four decades since the discovery of J/w in 1974. 
During this era, many theoretical studies have been carried out in the quark model framework to produce 
the spectrum of heavy quarkonium systems [4, 5]. An important manifestation of studying heavy QQ 
is that the spectrum of these states can be explained by using nonrelativistic formalism. For instance, 
the Cornell potential model [6, 7], which incorporates a spin-independent color Coulomb plus scalar 
linear confined potential, was hugely successful at describing the spectrum of the charmonium systems. 
The Cornell spin-independent potential is an approximate heavy quark spin symmetry (HQSS) within cé 
systems. Deviations with experiments can be observed in such potential models due to HQSS breaking 
effects [8, 9]. One possible source of breaking HQSS is the spin-dependent potential which introduces 
relativistic corrections to the Cornell potential model. The widely used relativized quark potential model, 
sometimes also referred as the Godfrey-Isgur (GI) model [10], is so far considered the best available option 
to reproduce the spectra of heavy quarkonium systems. 

In heavy QQ systems, hadronic transitions serve as a crucial probe of their internal structures and 
help to establish the understanding of light quark coupling with a heavy degree of freedom. In QCD, the 
well-established formalism for hadronic transitions is multipole expansion (ME) [11, 12, 13, 14], which 
assumes that these transitions take place due to the intermediate process of gluon emission. These gluons 
are supposed to be soft, having wavelengths much larger than the size of a heavy quarkonium system. 
These soft gluons further couple to 7(s) and 7 to complete such kinds of hadronic transitions. 

Development of heavy meson chiral Lagrangians (HMCL) [15] is the foremost simplification to QCDME. 
HMCL serve as an effective field theory (EFT) to QCDME in a soft exchange where the gluonic exchanges 
are predominantly of limited momenta. With the assumptions that (i) the heavy QQ involved in the pro- 
cess is well separated to consider it in a stringlike picture and (ii) the momentum of the emitted light 
meson is not too large, the HMCL are successful at reproducing the hadronic transitions among lower 
charmonia [16, 17]. The experimental status of the spectrum of higher vector charmonium(like) states is 
very rich now and several precise measurements have been recorded for their hadronic transitions [18]. 
To describe the observed transitions of higher cE systems there is a potential need for a theoretical model 
which can predict the transitions in the high momentum regime and help to identify the missing higher 
states through their hidden-flavor decays. We try to fulfill this need by modeling hadronic transitions of 
higher vector charmonia. Our proposed model is away from all the assumptions [(i) and (ii)] of HMCL 
and useful to predict the transitions involving much large momenta. 

Another possibility is that the transition between two S waves, S to P or D to S wave charmonia with 
the emission of 7 (7) might occur through intermediate open-charm contributions. Heavy quarkonium 
states can couple to intermediate heavy mesons through the creation of a light quark-antiquark pair. 
The formalism which incorporate intermediate heavy mesons within hadrons is sometimes referred to as 
coupled-channel effects. For instance, using the 3P) quark pair creation mechanism [19], intermediate 
meson loop contributions are found to be essential to explain the suppression of dielectric decay widths 
of higher bottomonium [20]. Coupled-channel effects have also been taken into account in the QCDME 


framework to study the hadronic transitions with two-pion emission for the charmonium system and found 
a good agreement with experimental measurements [21]. In this paper, we neglect the coupled-channel 
effects for simplicity, which can be included in the future in the unquenched quark model. 

To investigate the intermediate charmed meson loop effects on (3686) + J/wn decay, nonrelativistic 
effective field theory (NREFT) formalism was constructed [9, 22, 23]. It is noted that if we go to much 
higher waves e.g., Y(nS) or y((n — 1)D) with n = 4,5,6..., the decay momentum is not so small, as it 
lies in the relativistic regime; hence, the NREFT formalism is not very suitable for studying hadronic 
transitions of higher charmonia. 

These indications bring out the fact that there is a need for a comprehensive model which is capable of 
producing hadronic transitions with the emission of light meson(s) for higher mass charmonium systems. 
We attempt to fill this gap by modeling the transitions UV > J/yn and V > h,(1P)n, where W refers 
to n°S, and (n — 1)°D, vector charmonia with (n = 2,3,4,...). We present our predictions for hadronic 
transitions of higher vector charmonia into J/w7n and h,(1P)n, which provide useful references to determine 
their properties in ongoing and forthcoming experiments. 

In hadron physics, the most widely used model to study open-flavor strong decays is 3P) or the quark 
pair creation (QPC) model. Within the framework of the ?Po model, quark pair creation is induced from 
QCD vacuum. Hence, the generated quark pair shares the quantum numbers of vacuum (JPC = 0*++); 
therefore, it is referred to as the 3P) pair creation mechanism. The traditional 3P) model has been widely 
used in hadron spectroscopy and decays [20, 24, 25, 26]. In the °P) model, the probability to generate 
qq pairs is independent of the distance of the generation point from the valence quarks. In this work, 
we introduce the pair creation triggered from the Nambu—Jona-Lasinio (NJL) four-point-like effective 
interaction (LNL). Since Cygzy is a mixture of scalar and pseudoscalar interactions, it raises the quantum 
numbers of the created qq pair as a mixture of °P) and !Sp. It should be noted that the dynamics of the 
creation of a 3P) vertex in the NJL framework is totally different from the conventional °P) mechanism. 

Recently, the interaction of L£vjz, has been used to study the mixing of Q(sss) baryons with its 
pentaquark partner states and it is found that the Cnyr leads to strong mixing between three-quark and 
five-quark sss © sssqq (where q = u,d,s), with J? = 3, It was reported that this expected mixing 
results in the decrease of the energy of the lowest state [27]. There is no hint of such kind of mixing 
within the conventional 3P) model (which only involves the scalar interaction). Hence these are charming 
motivations to consider this interaction to study the hadronic decays of higher quarkonia. 

The paper is organized as follows. In Sec. 2, we review the development of the NJL model and its 
application in hadron spectroscopy and decays, where we deduce the effective Lagrangian for hadronic 
transitions with the emission of light meson(s). We conclude this section with an overview of well- 
established S — D mixing formalism. Section 3 is devoted to discussing the results for r(Y% > J/yn), 
I(v + h.(1P)n), and estimates of the branching fraction B(Y — Jọyn) for Y(4360), Y (4390), and 
Y (4660). Finally, we give a short summary in Sec. 4. 


2 Theoretical Framework 


2.1 NJL Motivated Effective Lagrangian 


Effective field theories are very useful when the dynamics of the system involves only a few relevant degrees 
of freedom instead of all. The NJL model is one of the best examples of such kinds of effective theories 
which have the capability to recover almost all of the features of the exact leading theory. Historically, 
Nambu and Jona-Lasinio modeled a scheme to explain the pions as nucleon-antinucleon bound states 
[28, 29]; afterwards, the scheme gained more appreciation for being used at a more microscopic level by 
changing the nucleon field into a quark field Y. NJL model only involves the quark degree of freedom, while 
the gluon degree of freedom is frozen in its point-like interaction vertex. NJL four-point-like interaction 


between quarks can be described by 


Ne 


Lrg = Egs XO [CPX YW)? + (PAi?) (1) 


2 a=0 


where Ne = 8 indicates the color degree of freedom; A*(a = 1,---,8) are Gell-Mann matrices in SU(3) 
flavor space with flavor singlet \° = 32, where Z is the unit matrix in the three-dimensional flavor 
space, Y represents the quark field and gs is the coupling strength. This four-point-like color interaction 
with only one free parameter (gs) has the capability to produce quite good results for the spectrum of 
low-lying and excited light mesons [30, 31]. Using the NJL motivated SU(2)@SU(2) chiral Lagrangian 
for the excited pions, p and w mesons, the strong decay widths for the V’ > PP, V! > VP, P! > VP 
transitions (V’ and P’ are the excited vector and pseudoscalar meson decaying into the vector V and 
pseudoscalar P meson, respectively) have been computed and found to be in good agreement with the 
experimental data [32, 33]. During the 1990s, attempts were made to extend this approach to study the 
radiative transitions and strong decays of the charmed mesons. It is noted that the qualitative estimates 
of the strong decay widths using this approach agree well with the experimental data [34]. 

In light of these phenomenological studies we propose another quark pair creation mechanism which 
is inspired by the NJL model. We model the coupling of the light scalar and pseudoscalar meson with the 
charm quark. The effective Lagrangian of our model contains both the scalar and pseudoscalar interactions 
as present in the NJL model. In the case when we link the light qq production with an (anti)quark line, 
the effective Lagrangian of our proposed model can be written as 


Lr= (by <o>+pivy <1 >), (2) 


where g is the overall coupling strength, ~ is the heavy quark field, and < o > and < 7 > are SU(3) 
singlet scalar and pseudoscalar meson, respectively. Since wy is the SU(3) singlet, the light sector should 
also be in a singlet. That is why the above Lagrangian does not have SU(3) flavor matrices as present 
in LynyjL. The color index can also be suppressed. The above Lagrangian £z allows the coupling of the 
(anti)quark line only to a scalar or isospin singlet pseudoscalar. The possible Feynman diagram for the 
process UV — J/w7 is shown in Fig. 2. 
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Figure 2: Quark level diagram of higher vector charmonia decaying into J/wn. 


In principle, 7 can also couple to the antiquark line, so one needs to calculate two diagrams. But going 
through the details, it becomes clear that both diagrams are equally contributing. 

Since the experimental data are available for UV > J/wn, it provides us with an opportunity to test 
this model. No experimental data are available for the emission of a light scalar meson, so we mainly 
focus on 7 emission transitions to check the validity of the model. It is quite possible that our proposed 
model can be extended to produce di-pion transitions. It would be through the intermediate production 
of a ø meson, which further decays into tt z~. But it will involve final state interactions (FSI) between 


mt and m~ as investigated in Ref. [46]. This would be interesting but leads to intensive work which we 
will consider in our future studies. 

To calculate the hadronic matrix elements, we prefer the mock hadron prescription [35] to express the 
initial and final meson wave functions. The mock hadron is defined as a collection of free quarks with the 
wave function of the bound quarks in a physical hadron normalized to the physical mass. The advantage 
of using this prescription is that it lets us calculate hadronic amplitudes as integrals over free-quark 
amplitudes. The initial mock meson wave function can be expressed as 


JA) = V2Ba $ (bm, $8.|Jamya) | poalla): B-P); (3) 
LS 


where E4 is the total energy of the meson; |Lm), |S.S), and | Jama) are the orbital angular momentum 
between the quark and antiquark, the total spin of the quark-antiquark pair, and the total angular 
momentum of the meson, respectively; (Lm, SS.|J4m,) is the Clebsch-Gordan coefficient; and ¢4 and 
x? are the spatial and spin wave functions of the initial meson |A), respectively. The relative momentum 
between the quark and antiquark, pj, is integrated over all values. We adopt the relativistic normalization 


. (A(p.)|A(p)) = 2E48° (p — p”). (4) 
j- For two-body A —> BC decay, we define the transition amplitude as 
O) (BO|H1|A) = 2r y 8E4EgEcô* (pi — pp) M™a™ Bc, (5) 


Considering the standard relativistic phase space, we define the decay width in the center-of-mass (CM) 
frame as 


EgE 
Tasso = 2rk A Y | Qg Mae], (6) 
MA MIp NIG 
= for any fixed mg}. Since the decay width is independent of the polarization of the initial state, we set 
© my, = Ja in the following calculations. Here k expresses the momentum of the outgoing mesons B or C, 
a which is given as 


vlna (me ~mo}*llm = (me + mo)? 


7 
— | (7 
C with Eg = \/m3, + k? and Ec = 4/mĉ + k?. The overlap of the wave functions of the initial meson |A) 


and the final mesons |B, C} can be expressed as 


M™a™ig™Ig = g | Emaoh: 5 tpPp)Mo, (8) 


where tp = mQ/(meg + mg) = 1/2, with mo = mg = Me, with me the charm quark mass. Mo is the 
free-quark amplitude which is discussed in detail in the Appendix A. To find the overlap of the wave 
functions, we use simple harmonic oscillator (SHO) wave functions, which can be written in momentum 
space as 


Wn,im(P) = Rn, 1(p)Vi" (p, 0, p), (9) 


where ny, l, and m represent the radial, orbital, and magnetic quantum numbers, respectively. Vi” (p, 0, y) = 
ply" (0, p) is the solid harmonic defined as a function of spherical harmonic. The radial wave function 
R,,1(p) is given as 

2n,! 


Tee er ras oP 6%, (10) 
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Rn,1(p) = 


i4i : : : 
where is an oscillatory parameter and Ln? (p?/ B>) is the associated Laguerre polynomial. Indeed, SHO 
wave functions serve as a coarse approximation to the true wave functions. However, qualitatively, SHO 
wave functions are similar to the realistic wave functions and useful for producing analytical results. 


2.2 S—D Mixing 


It is predicted that the J”? = 177 charmonia near or above the open-charm threshold are an admixture 
of S and D waves [36, 37, 38]. An S wave dominant state has a D wave component in its wave function 
and vice versa. A well-established formalism of this S — D mixing is based on reproducing the dielectric 
decay widths to deduce the mixing angle. If we neglect the open-charm contributions due to coupling to 
corresponding decay channels, under the assumption that the n?.$; state only mixes with (n — 1)? D1, S 
and D wave dominant states can be expressed as 


phys = cos O[n* $1) + sin 6|(n — 1)°D1) (11) 


Ups = —sin6|n°$,) + cos | (n — 1)? D1) (12) 


Here {phys and Wiss represent the S wave and D wave dominant state, respectively. The relative sign 
between |n°$,) and |(n — 1)3 D1) is just a matter of convention. One may follow the other convention as 
in Ref. [37], but the effect of relative sign can be compensated by swapping 0 + —0. A rough estimate of 
the S — D mixing angle can be made by computing the ratio of the dielectric decay widths [36, 37]. This 
has been done in the quark model framework by computing the wave functions using a potential model 
and then tuning the mixing angle to reproduce the dielectric decay widths. Considering (3770) as the 
1D dominant state with small 25 component, there exist two sets of possible ranges for the values of the 
mixing angle: 6 ~ —10° ~ —13° and 6 ~ +26° ~ +30° [36, 37, 39]. 

There also exist a couple of quark-model-based phenomenological studies in favor of large S — D 
mixing [38, 40]. A large mixing angle such as 0 = 34° is used [38] to produce almost the same dielectric 
decay widths of 7(4040) and 7(4160), which is consistent with experimental measurements [18]. The 
experimental fact that the n°.S; dominant states have relatively small dielectric decay widths while n? Dı 
dominant states have rather large widths can also be well described by considering the large S — D mixing 
as the underlying mechanism [38]. 


3 Results and Discussions 


3.1 T(Y¥ > Jn) 


To compute the decay widths for the process V —> J/wn, it would be better to analyze the dependence of 
the wave function on the oscillatory parameter ĝ in the case of SHO wave functions. For light qq systems, 
the best value of 8 from spectroscopy and decays is 0.379! [24, 41]. But for heavy QQ states it is bit 
larger and also has a range 8 = 0.4 ~ 0.6 in the literature [25, 42, 43, 44, 45]. The parameter 8 relates to 
the size of quark-antiquark bound state. Since the size of heavy QQ is smaller than the light qq system, 
Bag ~ Aqcp for q = u,d,s and Bgg ~ mau for Q = c,b, where mg is the mass and v is the velocity 
of the heavy quark. To choose the same ( for initial and final mesons is quite reasonable [46], although 
there exist some predictions that each meson has its own effective 3 [47]. Quark model studies show that 
the effective 8 for higher cé multiplets is smaller than the corresponding lower ones. It is due to the fact 
that the excited states have large spatial extensions. For instance, for 1P multiplets, 6 = 0.514 and for 
2P multiplets, it is 0.435 [42]. This is an indication that for higher charmonium states, the favorable 
value of 8 should be around 0.4. The value 8 = 0.44 has also been used to incorporate the spin counting 
predictions of open-flavor strong decays of higher S wave cé states under the °P) framework [44, 45]. We 
tune the parameter 8 along with the coupling strength g of the model. The coupling constant g is fitted 
by choosing 8 = 0.4 for all involved mesons, which agree with recent similar studies [43]. In principle, 
one can choose different values of 8 for initial and final mesons to be more accurate, but for simplicity we 
choose the same £ for all charmonium states. 


16 is in units of GeV; however, for the sake of simplicity, only its numerical values are written. 


Me = 1.50 GeV B=040 GeV g=0.80 |6| = 13° 


Table 1: The parameters used in our calculation. Due to the implicit treatment of color and flavor degrees 
of freedom, these factors do not show up in our calculation. 
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Figure 3: (Color online) 6 dependence of decay widths of the first few W’s into J/yn. The dashed black 
curve represents the pure S and D wave decay width; the solid blue and dotted red curves represent the 
decay width with small and large mixing angles, respectively. It is a coincidence that the decay width 
for W(3686) > J/wW7n is roughly the same for 6 = 13° and 0 = 34°, which causes an exact overlap of the 
curves. It is due to the definition of the S — D mixing mechanism. To confirm this argument, a rough 
estimate for a specific value of 8 can be made from Fig. 4, which contains the decay width dependence 
on the mixing angle. 
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Figure 4: (Color online) 6 dependence of the decay widths of experimentally well-established higher vector 
cé states decaying into J/y7. The dashed blue line represents the pure S or D wave, the dotted-dashed 
orange line represents 0 = 90°, and the solid red curve represents the decay width as a function of the 
mixing angle. The gray band indicates the experimental values of [(Y — J/pn) with statistical errors 
listed in Table 2. 


We explicitly show the 6 dependence of the decay width in Fig. 3 for ~(3686), (3770), (4040), 
(4160), and w(4415), to clarify the possible acceptable range of the 8. It is worthy noting that the decay 
width does not change drastically around 8 = 0.4. It is clear from Fig. 3 that one can choose any value 
of 8 within the safe region, i.e., 0.3 ~ 0.6. Our preferred value 8 = 0.4, lies in the safe region and, hence, 
is perfectly adequate. 


State  n?SHIL; ~— Trotai [18] Bip > J/yn) [18] Toos TSS [18] 


y>JI/ pn 
(3686) 23 Sy 0.296 + 0.008 (3.36 + 0.05)% 0.010 0.010 + 0.001 
(3770 1° Dı 27.2 1.0 94x 1074 0.025 0.025 + 0.011 
(4040 3°54 80+ 10 5.2+0.7 x 107° 0.347 0.416 + 0.076 
w(4160 23D, 70+ 10 < 8x 1073 0.204 < 0.560 + 0.080 
(4415 43S) 62 + 20 <6x 107% 0.425 < 0.372 + 0.120 


Table 2: All the widths are in units of MeV. For expressing the quantum numbers, we use the spectroscopic 
notation n25+!'L 7, where n = n, + l; n, is the radial quantum number; and S, L, and J represent the 
spin, orbital, and total angular momentum of charmonia, respectively. 


The parameters used in our calculations are listed in Table 1. Table 2 shows the fitted results with 
the choice of best-fit values of the parameters of Table 1. We get quite impressive agreement with the 
experimental data. Although there exist only upper limits for the ~(4160) > J/yn and 7(4415) > J/yn 
decay processes, our computed decay widths for the former decay process lie within this limit, while for 
the latter process our predicted width is slighter larger than the central value. It is worthy noting that 
the experimental value of [(w(4415) + J/wn) has large statistical errors. Considering this error range, 
our prediction in this case still lies within the upper limit. 

Error estimation in the theoretical model is still an open question. It became an important debate 
among theoretical constructs in the last few years. In general, there are no surefire prescriptions for 
assigning error bars to theoretical models. In the case of parameter dependence, by doing numerical 
analysis one can confine the model space to a physically reasonable domain. Within this domain, there is 
a range of reasonable parametrizations that can be considered as delivering a decent fit [48]. Using this 
prescription we scan the parameters of our model and find the physically reasonable range of 6 = 0.3 ~ 0.6 
GeV and |0| = 10° ~ 13°. For giving an idea of the uncertainties arising from model parameters, we 
explicitly give the decay widths in Table 3 by varying @ and 0 in the described physical range. 


Dyan 
State | 86=0.35 B=040 8=0.45 B=0.50]| |6|=10° jo] = 11° 0 =12 Pts | TSE y, [18] 
4p(3686) 0.017 0.010 0.006 0.004 9.782 x 107? 9.865 x 107? 9.942 x 107° 0.010 0.010 + 0.001 
w(3770) 0.042 0.025 0.016 0.011 0.041 0.035 0.030 0.025 0.025 + 0.011 
(4040) 0.906 0.347 0.146 0.066 0.336 0.340 0.344 0.347 0.416 + 0.076 
w(4160) 0.505 0.204 0.089 0.041 0.290 0.260 0.231 0.204 < 0.560 + 0.080 
(4415) | 1.647 0.425 0.123 0.039 0.411 0.416 0.421 0.425 | < 0.372 + 0.120 


Table 3: All the widths are in units of MeV and rounded to 0.001 MeV. While varying any parameter 
others are fixed and given in Table 1. 
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Figure 5: (Color online) Decay width dependence of Y — J/w7 on initial mass. Dashed black and red 
solid curves represent the decay width of S and D wave dominate states, respectively. 


The dependence of the decay widths on the S — D mixing angle 0 is very crucial to understand the 
behavior of this hadronic transition. We show the 6 dependence of r(Y — J/yn) in Fig. 4 along with a 
band gap, which actually represents the decay width range when we consider the corresponding w state 
as pure S and D wave states. Due to sinusoidal behavior, small mixing angles are adequate, while large 
mixing angles may ruin the predictions. Although Fig. 4 contains the predictions with a specific value of 
B, still it can be seen that the decay widths of the process 7)(3686) > J/wWn at 6 = 13° and 0 = 34° are 
exactly the same, which causes an exact overlap of the two curves in Fig. 3. 

We also give the estimates of the decay widths of ~(3°D1), (4°51), Y(43 D1), ¥(5°S1), Y(53 D1), and 
w(63.91) states decaying to J/wn in Fig. 5, which provide useful information to search and understand the 
missing higher vector states. Due to the fact that these higher states are poorly understood experimentally, 
we are not able to predict the widths exactly. To plot the decay width as a function of the mass of the 
corresponding higher cc vector state, we consider the mass range based on serval quark model predictions 


10 


of the mass spectra listed in Table 4. 


State JPO Screened [49] BGS-NR [50] BGS-Rel. [51] CQM [52] RS [38] 


w(4°S1) 177 4273 4406 ? 4356 4389 4420 
0(3°Di) 177 4317 4455 4470 4426 4470 
w(5°S,;) 177 4463 4704 4661 4614 4655 
(4D) 177 =- 4770 4735 4641 4700 
(6S1) 177 4608 4977 4912 4791 4815 
w(5°Di) 177 = — 4976 4810 = 


Table 4: Quark model predictions of mass spectra for higher vector charmonia. Reference [49] predicts 
mass in a nonrelativistic screened potential model, which incorporates the color screening effects due 
to the creation of a light gq pair within the heavy QQ. Reference [50] contains the predictions of the 
nonrelativistic effective QQ potential as calculated by Barnes-Godfrey-Swanson (BGS) [25]. Predictions 
of BGS potential with relativistic corrections to mass (calculated by using leading-order perturbation 
theory) are taken from Ref. [51]. Predications of the constituent quark model (CQM) [52] and Salpeter 
equation with relativistic string (RS) Hamiltonian [38] are also listed. “—” indicates that the prediction 
is not available. 


We include estimates of pure higher S and D wave states along with the predictions with small 
S — D mixing. In all considered cases, the decay width of the D wave state is smaller than that of the 
corresponding S wave one. Given that the D wave interferes destructively with the S wave, its decay 
width is suppressed significantly with the use of a small mixing angle. After a particular mass value, the 
decay width of both S and D wave states becomes very sensitive to the initial mass. For example in 
the w(3D/4S) case, the decay width rises exceptionally after 4.5 GeV. This critical mass value for other 
higher states can be observed from Fig. 5. With reference to the observed order of the 7 emission rate, 
these critical mass values might provide the upper limits on the mass of the corresponding higher vector 
charmonium. 


3.2 Predictions for ¢(nS/(n —1)D) > h.AP)n with (n = 3, 4,5) 


The hidden-charm n decay of JPC = 177 higher charmonia into the lowest P wave cé meson, i.e., 
h-(1!P,), is also possible. The threshold for this decay process is 4073 MeV. So, the first vector state 
which can decay into h-(1P)n is w(4160). There exists an experimental evidence for ete~ > h-(1P)n 
around 4170 MeV mass by the CLEO Collaboration [53]. Their reported measurement on the branching 
fraction B(7)(4160) > h.(1P)n) is < 2x 1073 with a 90% confidence level. The HQSS violating transition 
Ww — h-(1P)n requires the spin flip to be significantly suppressed relative to the corresponding heavy quark 
spin-conserving transitions like V > J/¢am [54]. The observed ratio r(¥ > J/yn)/T(¥ > J/yrr) is 
fully consistent with the earlier theoretical predictions [55]. It has been argued in [23, 56] that the 
coupled-channel effects due to intermediate charmed mesons for these transitions are quite small. 

Among well-known higher vector charmonia, only 7(4160) and (4415) have enough phase space to 
decay into h,(1P)n. Table 5 contains our predictions for these states. We predict the decay width of 
(4415) > h-(1P)n with the same order of magnitude as w(4160) to the similar final state: 


T(Y(4160) > he(1P)n) 
T(Y(4160) > J/Yn) 


= 7.887 x 107”, (13) 


?This value is taken from Ref. [25]. 
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T(Y(4415) > he(1P)n) 
D(p(4415) > J/Yn) 


It is not easy to give the estimates of the decay width for y(3D), y(4D), and (55) or higher ones 
because these states have not been experimentally well established up to now and hence, their masses are 
unknown. We give the initial mass dependence of the decay width of the V + he(1P)n transition of these 
higher vector states in Fig. 6, both for the pure S and D wave and for the standard S — D mixing case. 


= 6.736 x 107°. (14) 


State n?S**L; Trotai [18] Bly > he(1P)n) [53] Pn apn Tohe P) 53] 
(4160) 2° Dı 7010 <2 x 107° 1.609 x 107? < 0.140 0.020 
(4415) 43S: 62 + 20 — 2.863 x 107? — 
Table 5: All the widths are in units of MeV. “—” indicates that the experimental data are not available. 


__ WGDIAS), 0=13° 


a 


4.3 44 45 46 47 48 4.9 4.3 4.4 4.5 4.6 4.7 4.8 4.9 


Mass(GeV) Mass(GeV) 
y(4D/5S), 0=0° y(4D/5S), 0=13° 


Mass(GeV) Mass(GeV) 


Figure 6: (Color online) Decay width dependence on the initial mass for a few W’s into h.(1P)n. Dashed 
black and solid red curves represent the decay width of S and D wave dominate states, respectively. 


For Y —> h-(1P)n, the decay width of the pure D wave is larger than that of the corresponding pure 
S wave. It is due to that the overlap of the wave function of the D wave with 1P is larger than the 
corresponding S wave. The D wave interferes destructively with the S wave; the decay width of D wave, 
in this case, is not suppressed significantly as we have observed for the UV — J/wn channel. For the 
Y — J/wn case, the reason why the decay width of the pure D wave is smaller is that the overlap of the 
D wave with 1S is smaller than the corresponding S wave state. The swapping of the curves can be seen 
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in Fig. 6 with the increase of the mixing angle. The constructive interference of the S wave with the D 
wave causes a significant increase in the decay width of the S wave dominate state. 


3.3 Y(4360), Y(4390) and Y (4660) Assignments 


Just after the experimental observation of Y (4360) and Y(4660) in the initial state radiation (ISR) 
ete” + ysr TT W(2S) process at Belle [57], many theoretical studies were carried out to incorporate 
these states into conventional and exotic cé spectra (for an overview, see the discussion in Sec. 4.8 of the 
recent review [58]). Among these configurations, there exist a couple of interpretations by considering 
Y (4360) and Y(4660) as canonical 33D, [50] and 5.9; charmonium [50, 52], respectively. Predictions 
were also made for the dielectric widths, E1 and M1 transitions, and open flavor strong decays. In the 
screened potential model, the state Y (4360) is also interpreted as %(33 D1), while Y (4660) was considered 
as w(6°S;) [49]. Assuming Y (4360) as 33D, and Y(4660) as 5°.S, dominate, their dielectric decay widths 
can be reproduced to get agreement with experimental data by introducing the large S — D mixing [38]. 

The predicted mass of 7)(3°D,) in the BGS and RS potential model (as shown in Table 4) is somewhat 
larger than the experimental mass of Y (4360). However, we notice that the mass predictions of various 
potential models for the higher cé states may differ by 10 ~ 100 MeV [4, 59]. Therefore, it is not irrational 
to treat Y (4360) as w(3D) dominant with a small component of ~(4S), exactly in the same way as for the 
lower D waves in Sec. 2.2. Mass predictions for (5°51) of various potential models do not differ much 
from Y (4660) except the one produced by the screened potential model. Hence, it is also not irrational 
to treat Y (4660) as w(5S) dominant having a small component of 7)(4D). 

Very recently, the BESIII Collaboration announced the observation of a new resonant structure, 
Y (4390) in the ete~ — ntn“ he process [60]. Its measured mass and total decay width are (4391.6 + 
6.341.0) MeV and (139.5+16.1+0.6), respectively. The width of Y (4390) seems to be broader. It might 
be possible that this higher mass region could be described either by a Y (4360) resonance, or by a phase- 
space background as has already been noticed in the case of Y (4290) [61]. Since the reported quantum 
numbers of this state are JPC = 177 and it lies in the mass region of 3D charmonium, hence, this new 
resonance may also be considered as a candidate for 7(3D). As an estimate, we give our predictions of 
its 7 transition branching fraction in Table 7 by assigning 3° D1. 

One may argue that the vector state Y (4260) can also be assigned as a 33 D1 or 4°; charmonium state. 
As listed in Table 4, in various potential models, %(33 D1) and ~(4°S) are predicted between 4.3 ~ 4.5 
GeV. Therefore, with reference to the observed mass of Y (4260), it lies below the (45) ~ y(3D) mass 
region. The value of its [.+,- is two to three orders of magnitude smaller than those of the well-established 
conventional vector charmonia [62]. It is nearly impossible to accommodate Y (4260) as a conventional 
charmonium state. Many theoretical constructs consider Y (4260) as an exotic charmonium [61], for 
example, as a DD,(2420) molecule with a binding energy of 29 MeV [63], or as a heavy hybrid meson 
with a gluonic excitation of about 1 GeV higher than the lowest ne and J/w [58, 62]. 

Despite the fact that these Y states do not decay into open-charm channels, it would be interesting 
to study their hidden-charm strong decays. By assuming Y (4360) and Y (4660) as ~(3°.D 1) and (5°S;) 
dominant states, respectively, we give our predictions for Y (4360) > J/wWn and Y (4660) > J/wn, which 
might be helpful to understand the properties of these vector states. Because only experimental upper 
limits [64] exists for the product of the branching fraction B(Y — J/wn) and T.+,.-(Y) for Y (4360) 
and Y (4660), we need to know the dielectric decay width of these states. There exist few theoretical 
predictions for [.+,- for these states as summarized in Table 6. 


13 


Initial state n?St1 L; Mass (MeV) Te+e- (keV) 


T(Y (4360) +ete~) 33D, 4361 0.87 [50 
4455 0.83 [50 
4426 3 0.33 [52 
4470 0.06 [38 
4470 0.63 + [38] 
T(Y (4660) — ete) — 5391 4664 1.34 [50 
4704 1.32 [50 
4614 0.57 [52 
4655 0.73 [38 
4655 0.39 + [38] 


Table 6: ete~ dileptonic partial decay widths for Y(4360) and Y (4660) states. 


As shown in Table 7, the only available experimental information is B(Y > J/wn) - Des for these 
Y states. Therefore, we need [.+,- to compute [(Y — J/yn). Table 6 contains different theoretical 
predictions for l e+e- of these Y states both for pure S and D waves, and with the large mixing (@ = 34°). 
For Y (4360), we take the average value of [.4,- = 0.523 keV from Table 6 by considering it as pure 
3D. The experimental upper limit for T (Y (4360) —> J/w7) is given in Table 7. To give a comparison 
with Ref. [38], we also compute the upper limit of r(Y (4360) > J/wn) with large S — D mixing, i.e., 
0 = 34°. In all three cases, pure 3D, small, and large mixing, our predictions are in agreement with the 
experimental measurements. We conclude that Y (4360) could be considered as a potential candidate for 
dominant 33D; charmonium state. 


Ty J/yn DYS, s/n 
State n®>*1L, Tiotal B(Y > J/yn) | @=0° 0=13°  0=34° 0=0° 6@=34° 
Y (4360) 3D 74+ 18 [18] ros [64] | 0.047 0.016 10x107? | <0.963 < 0.799 
Y(4390) 3Dı 139.5 + 16.1 [60] z 0.083 0.028 1.6 x 107° = -— 
Y(4660) 53S; 48+15 [18] p24 [64] | 0.057 0.070 0.077 | <0.046 < 0.116 
Table 7: Predictions for r(Y — J/wn) for the Y (4360), Y (4390), and Y (4660) states. “—” indicates that 


the experimental data are not available. All the widths are in units of MeV. 


For Y (4660), we also include the predictions for pure 55 and the mixed case. For the pure 5°9; state, 
we take the average value of [.+,- = 0.99 keV from Table 6 and list the experimental upper limit along 


?Reference [52] assigned 3?D, to (4415) instead of Y (4360), while our concern here is just having a comparison among 
predictions of the dielectric decay width of 33D, with different input masses. 
“This prediction contains the admixture of 4°.5; for Y(4360) and 4°.D, for Y(4660) with mixing angle 0 = 34°. 
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with our prediction for I (Y (4660) —> J/w7) in Table 7. In the large S — D mixing case (6 = 34°), the 
dielectric decay width = 0.39 keV [38] allows us to give an upper limit on r(Y (4660) > J/yn). Our 
predicted value, in this case, is within this upper limit as shown in Table 7. In all cases, our predictions 
agree with the experimental data. Hence, our results are consistent with the experimental data and the 
state Y (4660) can be considered as 7(5°S,) dominant with a small y(43 D1) component. 

In the case of Y (4390), for the sake of completeness, we give our predictions of its hidden charm 7 
decay, with and without S — D mixing. To identify this state, measurements on its hadronic branching 
fraction are required. We think that these estimates might be useful to clarify the picture of these vector 
states and give some references to search for the missing higher S and D wave vector charmonia. 


4 Summary 


A model to create a light meson for heavy quarkonium transition is proposed. This model is used to 
study the decays of higher vector charmonia into J/yn and h-(1P)n. Computed decay widths are in 
excellent agreement with experimental data. The ratio [(V > h.(1P)n)/T(Y — J/wn) is predicted for 
(4160) and (4415). The initial state’s mass dependence of r(Y — he(1P)n, J/wWn) for higher vector 
charmonium is given. We suggest that the ongoing (Belle and BESIII) and forthcoming (PANDA and 
BelleII) experiments look for suggested unobserved decay channels. We also give the estimates of 7 
transition branching fractions for Y (4360), Y (4390), and Y (4660) by assuming them as cc bound states 
with quantum numbers 3? D1, 33 D1, and 53.51, respectively. Our predictions reflect that the state Y (4360) 
can be considered as a potential candidate for the 3°D, charmonium state. Assuming Y (4660) to be 
5°91, the predictions are consistent with the experimental upper limit. For a broader Y(4390) state, the 
update on its hadronic branching fraction from BESIII is eagerly awaited. We hope that our predictions 
might provide useful references to determine the properties of higher charmonium states in ongoing and 
forthcoming experiments. 
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Note added in proof 


After the revision of this manuscript, BESIII published an evidence of ete~ — nh, at center-of-mass 
energy ys = 4.358 GeV [65]. Along with their earlier measurement [66], this evidence will help to delve 
into the Y (4360) through its HQSS violating hadronic transitions. 
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A Free-Quark Amplitude 


To evaluate the matrix element of A — BC decay, we need spin matrix elements. At the quark level, 
these matrices involve the matrix elements of the Dirac bilinear (with T = iy? and J in our case) and 
Pauli matrix elements. Matrix elements of the Dirac bilinear in the nonrelativistic limit can be expressed 
as 


lim Gy Bsa m (15) 
1m Jol = p 
gota te ) t (s!lals)-(@- 7) T= i 

—dggr T = I 
lim ssl vyz = F 2s 16 
go 7 e -(-F) T=iy’. 06) 


We have already shown that it is easy to handle the wave function overlap integration in a Cartesian 
basis; therefore, we express the elements of Pauli spinors in terms of Cartesian basis vectors as 


(L |] t) = (ê +ib) 
(t ElL) = (@- 19) (17) 
(t lal t) =—-( lal) = 2 


For antiquark case, these relations are 


(lait) = -@-@) 
(ail) = -(@ +48) (18) 
lait) = -Ual = -2 


From Fig. 2, one can get the following relation of momentum conservation by considering the initial state 
in the center-of-mass reference frame: 


P2 = —P1 
pi =pi—k (19) 
Pp = —k. 


The meson’s space wave function can be written as 


TER (20) 
P = m 

OB{ UB ay = 7a = bp(pi — aie BD) 

where the m;i = (A,B) are the reduced masses of the constituent quarks of the mesons A and B, 


respectively. It should be noted that all the momenta should and can be expressed by the integration 
variable and the momentum of the B meson, viz. pı and Pg in this work. The free-quark amplitude for 
Fig. 2 is 

iMo = [tgs (iy Jgs] [Ugl vgs’), (21) 


i 


iMo = ;—(1'|ë]1) - (P1 — P1) (2165512), (22) 


2Me 


where me is the mass of the charm quark. Collecting all the pieces together, the full amplitude becomes 
a SEEE 3 5 ae 
Mamae = grt | Ppa PARP- BPa Fi — Pr) Qldsel2"). 23) 
© 


This is the nonrelativistic approximation of Eq. (8). The momentum Pg can be computed by using 
Eq. (7). 
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